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We present results of wavepacket simulations of scattering of an oriented methane molecule from 
a fiat surface including all nine internal vibrations. At a translational energy up to 96 kj/mol we 
find that the scattering is almost completely elastic. Vibrational excitations when the molecule hits 
the surface and the corresponding deformation depend on generic features of the potential energy 
surface. In particular, our simulation indicate that for methane to dissociate the interaction of the 
molecule with the surface should lead to an elongated equilibrium C-H bond length close to the 
^ , surface. 
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Oh 1 I. INTRODUCTION 

g: 

The dissociation of methane on transition metals is an important reaction in catalysis. It is the rate limiting 
step in steam reforming to produce syngas.El It is also prototypical for C-H activation in other processes. A large 
numbe* p£ molecular beam experiments in which the, dissociation energy was measured as a function of translational 



energymla as well as various bulb gas experimentsli^cil have already been done on this system. These experiments 
have contributed much to our understanding of the mechanism of the dissociation. However, our knowledge is still 
r/3 , small compared to what we know about a well-studied reaction as H2 dissociation. 

Some years agp-,we performed a multi-configurational time-dependent Hartree (MCTDH) study of CH4 dissoci- 
ation on NiLLLl^?3 with a potential energy surface (PES) based on our earlier density functional theory (DFT) 
calculations .olj We included the distance of the methane molecule to the surface, a C-H distance, and the orienta- 
tion of methane as coordinates. OtherpWaizEpacket simulations focused on (some of) these coordinates in combination 
with lattice motion on several metals.ETO Othyc-.theoretical studies on this system used DFT calculationsEjLD too 
and classical stochastic trajectory simulations EaEJ 

None of the wavepacket simulations published so far have looked at the role of the internal vibrations of methane. 
It was observed experimentally that vibrationally hot CH4 dissociates more readily than cold CH^wj-jjhptlie energy in 
the internal vibrations being about as effective as the translational energy in inducing dissociation.B0Q&tL3 However, a 
more detailed assessment of the importance of the internal vibrations could not be made, because of the large number 
7— I I of internal vibrations. A recent DFT calculation also showed that the transition state for CH4 dissociation on Ni(lll) 
involves considerable internal excitation of the moleculeEj In this paper we report on wavepacket simulations that 
we have done to determine which and to what extent internal vibrations are important for the dissociation. We are 
not able yet to simulate the dissociation including all internal vibrations. Instead we have simulated the scattering of 
methane, for which all internal vibrations can be included. By looking at vibrational excitations and the deformation 
of the molecule when it hits the surface we can derive information that is relevant for the dissociation. We have used 
>^ PESs that have been developed with Ni(lll) in mind, but our results should hold for other surfaces as well. 

We have again used the MCTDH methojifar the wavepacket simulations, because it can deal with a large number 
of degrees of freedom-and,with large grids. EJtlJ This method has been applied successfully to gas phase reactions and 
reactions at surfaces HaEia 

The rest of this paper is organized as follows. We start with a brief description of the MCTDH method. Then the 
various PESs that we have used are derived. (A harmonic intramolecular PES is adapted to include anharmonicities 
in the C-H distance, the decrease of the C-H bond energy due to interactions with the surface, and the increase of 
the C-H bond length also due to interactions with the surface.) The results of the simulations are presented and 
discussed next. We focus on excitation probabilities and deformation of the molecule when it hits the surface. The 
implications for the dissociation are discussed separately. We end with a summary and some general conclusions. 



II. COMPUTATIONAL DETAILS 
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A. The MCTDH Method 



We give here a short overview of the MCTDH approximation for completeness. More details can be found in 
Refs. |40| and 41. The exact wave- function of a D-dimensional system, is approximated by an expression of the form 

*MCT DH (gi,...,to;t) = J2 a 'n...n D (t)^(i^ t )---^i D o > (iD;t). (1) 

n 1 ...n D 

From this expression, it is possible to obtain the equations of motion for the one-dimensional functions ipn](qi]t) and 
for the correlation coefficients c ni ... nD (t). Without loss of generality we can choose the V , "i s to be natural single- 
particle states.cil They insure that we obtain the best approximation to the exact $f(qi, ■ ■ ■ , qr>\ t) for a fixed number 
of configurations; i.e., they minimize the expression (A|A) where 

A = * - *mctdh, (2) 

and where \& is the exact wave-function. The natural single-particle functions are eigenstates of the reduced density 
operators. 

Pjilj'l'j) = J d< li ■ ■ ■ dqj-idq J+1 . . . dq D *(<7i . . . ^•_ 1 g J g J+1 . . . q D )^*{qi ■ ■ ■ Ij-x^qj+i ■ ■ ■ to)- (3) 

The equations of motions for the natural single-particle states can be obtained by differentiation of the eigenvalue 
equation: 



dq' j p i {q j ,q'M%Xq' j ) = v$^{q j ). (4) 



This gives us 



dt 3 33 ^ 3 3 3 m^k]) t! 3 



^ril (5) 



where 



and 



mm -mm) ' 
= e e - ni ... nD ^...^^--m D \ (7) 

ni ...nj-i rij+i ...n D 

with the Hamiltonian H given by 

D 

The equations of motion for the coefficients 

Cn 1 ...nn=(^--^ i n D ) m (9) 

are again obtained by differentiation. 

D 

dt 



%h^-C ni ... nD = (< ) ...C ) |F|*)-Ecni...n 3 --im 3 -n j + l ...n ^Jn,. (10) 



Equationp-^) and ( |Toj ) are a particular form of the more general equations obtained from a time-dependent variational 
principleJUJ They conserve the norm of the wave-function and the mean energy of a time-independent Hamiltonian. 



2 



The resulting system of first-order differential equations, has to be solved with a generaLpurpose integrator. We 
used the variable-order variable-step Adams method, as implemented in the NAG libraryO This method gave good 
convergence for all described simulations. For the hardest simulation the method needed 195027 calculations of the 
time derivatives for 35000 atomic units simulation time, which took around two days CPU-time on an 166 MHz. 
SUN Ultra SPARC. The singularities in Eqs. (|5|), (||), and ( |Tc| ) have been treated numerically by the regularization 
procedure described in Ref. |4l| . 

The natural single-particle states have some small advantages over other possible choices of single-particle states. 
The most important one is that one can directly see from (ipn^ I'ipn^) how well $mctdh approximates the exact wave- 
function. How much a natural single-particle functions contributes to the wave-function is given by the eigenvalue Vn 
of the reduced density matrix. In an approximate MCTDH simulation {^n^\4>n^) is an approximation for this exact 
eigenvalue. The natural single-particle functions are also convenient for interpreting the results of a simulations. 

B. The Potential Energy Surfaces 

The PESs we used can all be written as 

Vtotal = Mntra + U sur f , (11) 

where Vmtra is the intramolecular PES and U S urf is the interaction with the surface. For the Vmtra we looked at four 
different types of PESs. 



1. A harmonic potential 

The first one is completely harmonic. We have used normal mode coordinates for the internal vibrations, because 
these are coupled only very weakly. In the harmonic approximation this coupling is even absent so that we can write 

Mntra &S 

1 10 

Mntra = Vharm = - £ k.Xf , (12) 
i=1 

the summation is over the internal vibrations, X^s are mass- weighted displacement coordinates and fcj are mass- 
weighted force constants. i— /see Table § for definitions and values); (X\ is the mass- weighted overall translation 
along the surface normal) .l3 The. force constants have been obtained by fitting them on the experimental vibrational 
frequencies of CH4 and CD4.E3E3 

We have assumed that the interaction with the surface is only through the hydrogen atoms that point towards the 
surface. We take the z-axis as the surface normal. In this case the surface PES is given by 

A ^ 

^ = ^E e " aZi ' ( 13 ) 

i — 1 

where Nh is the number of hydrogens that points towards the surface, a=1.0726 atomic units and ^4=6.4127 Hartree. 
These parameters are chosen to-give the same repulsion as the PES that has been used in an MCTDH wavepacket 
simulation of CH4 dissociation.^ 

If we write V sm { in terms of normal mode coordinates, then we obtain for one hydrogen pointing towards the surface 

V snxi = Ae- aiXl e-~ a2X2 e~ a3X3 e- aiX \ (14) 

where A as above, and a's as given in Table |n]. X2, A3 and A4 correspond all to a\ modes of the C^ v symmetry. 
There is no coupling between the modes A5 to Xiq in the V SUI f part of the PES, which are all e modes of the C3,, 
symmetry. Fig. |](a) shows a contour plot of the cross-section of the total harmonic PES with one hydrogen pointing 
towards the surface in the translational mode Ai and the v% asymmetrical stretch mode X4. 
For two hydrogens we obtain 

V BUI f = Ae~ aiXl e- a2X2 e- a3X3 e- aiX4 e- a5X5 (15) 

1 r 
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with A again as above, a's and /3's as given in Table ||. X2, A3, X4 and X§ correspond all to a\ modes of C% v . X7, 
Xg, Xg and Xiq correspond to b\ and 62 modes of C21,. Xq corresponds to the 02 mode of C2„ and has no coupling 
with the other modes in V SUI { . 
For three hydrogens we obtain 



V snlf = Ae^ aiXl e- a2X2 e' a3X3 e~ aiXi 
1 



(16) 



_g/3i^5 g — fcX 6 g— fe X 7 e 0iX & e l3sX 9 g—PeXto 



_ e - 2ft X 6 g2/3 3 X 7 g - 2/3 5 X 9 



X2, A3 and A4 corresponds to a± modes in the C3„ symmetry. Because these last six coordinates correspond to 
degenerate e modes of the C3„ symmetry, the (3 parameters are not unique. 



2. An anharmonic intramolecular potential 



Even though we do not try to describe the dissociation of methane in this paper, we do want to determine which 
internal vibration might be important for this dissociation. The PES should at least allow the molecule to partially 
distort as when dissociating. The harmonic PES does not do this. A number of changes have therefor been made. 
The first is that we have describe the C-H bond by a Morse PES. 



Morse 



d, 
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i=l 



(17) 



where D e — 0.1828 Hartree (the dissociation energy of methane in the gas-phase) and Ar^ the change in bond length 
from the equilibrium distance. 7 was calculated by equating the second derivatives along one bond of the harmonic 
and the Morse PES. If we transform Eq. ([l7]) back into normal mode coordinates, we obtain 



D, 
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g7i2^2 g7i3^3 g7i4-f 4 g7i7^7 g7i8^8 g7i9^9 g7», 10-^10 



(18) 



with D e as above. 7's are given in Tables III and IV. Note that, although we have only changed the PES of the bond 
lengths, the V4 umbrella modes are also affected. This is because these modes are not only bending, but also contain 
some changes of bond length. 

The new intramolecular PES now becomes 



Vintra — Vh ail n + Vuorsc ~ V COII , (19) 

where Vharm is as given in Eq. ([IJ) and V corl is the quadratic part of T4iorsc, which is already in Vharm- Fig. |l|(b) shows 
a contour plot of the cross-section of this total anharmonic PES with one hydrogen pointing towards the surface in 
the translational mode X\ and the 1/3 asymmetrical stretch mode X4. (We would like to point out that there are, 
of course, various anharmonic PESs for methane in the literature. There are two reascpSr-«fhy we haven't use them. 
First, these PESs are not in appropriate form to use then with the MCTDH methodxScil Second, these PESs are 
generally quite complicated. We prefer to keep it as simple as possible, because at this moment we're only interested 
in qualitative effects.) 



3. Intramolecular potential with weakening C-H bonds 

When the methane molecule approach the surface the overlap of substrate orbitals and anti-bonding orbitals of the 
molecule weakens the C-H bonds. We want to include this effect for the C-H bonds of the hydrogens pointing towards 
the surface. We have redefined the VMorse given in Eq. (|l^) and also replace it in Eq. ([H]). A sigmoidal function is 
used to switch from the gas phase C-H bond to a bond close to the surface. We have used the following, somewhat 
arbitrary, approximations, (i) The point of inflection should be at a reasonable distance from the surface. It is set to 
the turnaround point for a rigid methane molecule with translation energy 93.2 kJ/mol plus twice the fall-off distance 
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of the interaction with the surface, (ii) The depth of the PES of the C-H bond is 480 kJ/mol in the gas phase, but 
only 93.2 kJ/moi-near the surface. The value 93.2 kJ/mol corresponds to the height of the activation barrier used in 
our dissociations (iii) The exponential factor is the same as for the interaction with the surface. 
If we transform to normal-mode coordinates for the particular orientations, we then obtain 



4 2 

l/ wcak = D e 22 Wi 1 - e 7 * 2X2 e 7l3X3 e 7l4X4 e 7l7X7 e 7l8X8 e 7l9X9 e 7 ^ loXl0 , (20) 



i=l 

where Wi = 1 for non-interacting bonds and 

Wi = 7" e y ^ (21) 

for the interacting bonds pointing towards the surface. a\ is as given in Table || 7's are given in Tables HIJ and |IV| , 
Q = 1.942 ■ 10 _1 and u> = 7.197. Fig. |l](c) shows a contour plot of the cross-section of this total anharmonic PES with 
weakening C-H bond for the orientation with one hydrogen pointing towards the surface in the translational mode 
Xi and the ^3 asymmetrical stretch mode X4. 



4. Intramolecular potential with elongation of the C-H bond 

A weakened bond generally has not only a reduced bond strength, but also an increased bond length. We include 
the effect of the elongation of the C-H bond length of the hydrogens that point towards the surface due to interactions 
with the surface. We have redefined the VMorse given in Eq. ( Jig ) and also replace it in Eq. (|l9|) for this type of PES. 
We have used the following approximations: (i) The transition state, as determined by Refs. |27| , p8| , has a C-H bond 
that is 0.54 A longer than normal. This elongation should occur at the turn around point for a rigid methane molecule 
with a translation energy of 93.2 kJ/mol. (ii) The exponential factor is again the same as for the interaction with the 
surface. 

If we transform to normal-mode coordinates for the particular orientations, then we obtain 

V shiit = D e ^ W - e 7«JC2 e 7«^8 e 7<*^4 e 7<7XT e 7i8X Be 7 <9 -Jf9 e 7<,io^io e xp[S ie ~ QlXl ] , (22) 
i=l 

where a\ is as given in Table [n], 7's are given in Tables |Hj and |Vl For orientation with one hydrogen towards the 
surface we obtain; Si = 2.942 • 10 2 and S 2 = S3 = S 4 = 0, with two hydrogens; Si = S 2 = and S3 = S 4 = 1.698 • 10 2 , 
and with three hydrogens; Si = and S2 = S3 = S4 = 2.942 • 10 2 . Fig. |l|(d) shows a contour plot of the cross-section 
of this total anharmonic PES with elongation of the C-H bond for the orientation with one hydrogen pointing towards 
the surface in the translational mode X\ and the 1/3 asymmetrical stretch mode X4, and Fig. shows a contour plot 
of a cross-section in the translational X\ and the v\ symmetrical stretch mode X 2 in the three different orientations. 

Finally we like to note that the monotonic behaviour of all PES types in the translational mode does not contradict 
with the fact that the dissociative adsorption of methane is activated. The activation barrier for dissociative adsorption 
in the translational mode is situated in a region with very high excitations of the stretch modes, which we do not 
reach in these simulations. 



C. Initial States 

All initial states in the simulations start with the vibrational ground state. The initial translational part is 
represented by a Gaussian wave-packet, 

V> (tr) (A^i) = (27ra 2 r 1/4 exp 

where a is the width of the wave-packet (we used a = 320.248 atomic units), Xq is the initial position (we used 
Xq = llcr, which is far enough from the surface to observe no repulsion) and Pi is the initial momentum. We 
performed simulations in the energy range of 32 to 96 kJ/mol. We here present only the results of 96 kJ/mol 
(equivalent to P± — 0.2704 atomic units), because they showed the most obvious excitation probabilities for VMorse- 



(*i - x y 

4cj 2 



iPxXx 



(23) 
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We used seven natural single-particle states, 512 grid points and a grid-length of 15a for the translational coordinate. 
With this grid- width we can perform simulation with a ttsnslational energy up to 144 kJ/mol. 

Gauss- Hermite discrete- variable representations (DVR)E2l were used to represent the wavepackets of the vibrational 
modes. We used for all simulations 5 DVR points for the v 2 bending modes and 8 DVR points for the umbrella, v$ 
asymmetrical stretch, and v\ symmetrical stretch mode for an numerical exact integration, except for the simulations 
with Vshiftj where we used 16 DVR points for the v\ symmetrical stretch mode, because of the change in the equilibrium 
position. 

We did the simulation with one hydrogen pointing towards the surface in eight dimensions, because the v 2 bending 
modes X 5 and X 6 do not couple with the other modes. We needed four natural single-particle states for modes X 2l 
X3 and X4, and just one for the others. So the number of configurations was 7 1 ■ 4 3 • l 4 = 448. The simulation with 
two hydrogens pointing towards the surface was performed in nine dimensions. One of the v 2 bending mode (X G ) 
does not couple with the other modes, but for the other mode X5 we needed four natural single-particle states. The 
number of configurations was 7 1 • 4 4 • l 4 = 1792, because we needed the same number of natural single-particle states 
as mentioned above for the other modes. We needed ten dimensions to perform the simulation with three hydrogens 
pointing towards the surface. We used here one natural single-particle state for the modes X5 to X%o and four natural 
single-particle states for X 2 to X4, which gave us 7 1 ■ 4 3 • l 6 = 448 configurations. 

III. RESULTS AND DISCUSSION 

We found that the scattering is predominantly elastic. The elastic scattering probability is larger than 0.99 for all 
orientations and PESs at a translational energy of 96 kJ/mol, except for the PES with V^hift and three hydrogens 
pointing towards the surface for which it is 0.956. This agrees with the observation that the translation-vibration 
coupling is generally small.Ej 

If we would have wanted to determine the role of the internal vibrations from the scattering probabilities, we would 
have to do quite accurate simulations. We have opted instead to look at the molecule when it hits the surface, which 
enables us to obtain good results with much less costly simulations. 

A. Excitation Probabilities 

The surface PES has Cz v (with one or three hydrogens towards the surface) or C 2v symmetry (with two hydrogens 
towards the surface). If we expand this PES in a Taylor-series of internal vibrations, we see that the linear terms 
contain only those vibrations that transform as a± in C3„, respectively, C 2v . These are therefore easier to excite than 
others; We did not find any appreciable excitation of e modes of C3,, and the a 2 , bx, and b 2 modes of C 2v . 

We will not present results of the simulations with the harmonic PES, because they give almost the same excitation 
probabilities as the PES with VMorse- The maximum excitation probabilities at the surface for the PES with VMorse 
are given in Table We have observed the highest excitations for this PES in the V4 umbrella and v 2 bending modes 
in the orientation with two hydrogens, and for the umbrella mode in the orientation with three hydrogens pointing 
towards the surface. The excitation probabilities for the v\ and 1/3 stretch modes at this orientations are a factor of 
magnitude lower. 

We have observed higher excitation probabilities for the v\ and ^3 stretch modes in the one hydrogen orientation. 
The excitation probability of the umbrella mode is here lower then for the other orientations, but still in the same 
order as the stretch modes for this orientation. This can be explained by the values of the a's of V^ ur f (see Table 
||) and the force constants of Vq larm (see Table |), because Vharm is approximately VMorse for the ground state. The 
force constants of Vintra for the v\ symmetrical stretch stretch and z/3 asymmetrical stretch are of the same order, but 
the a parameter in V^ ur f in the orientation with one hydrogen pointing towards the surface for the u 3 asymmetrical 
stretch is around twice as large as for the v\ symmetrical stretch, which explains why 1/3 is more excited as v\. The 
surface repulsion on the 1/4 umbrella mode is even three times lower then on the v\ symmetrical stretch mode, but the 
force constant is also much lower. It results in a little more excitation of the umbrella then in the v\ symmetrical 
stretch. 

For the orientation with two hydrogens pointing towards the surface the repulsion on the vibrational modes is for 
all modes in the same order (v 2 bending a little higher), so the difference in excitation probabilities correlate here 
primarily with the force constants. The force constants of the v 2 bending and ^4 umbrella modes are of the same 
order, as are those of v\ and ^3. The stretch force constants are higher, however, so that the excitation probabilities 
are lower. The repulsion on the z/4 umbrella mode is the largest in the orientation with three hydrogens pointing 
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towards the surface. The force constants is lower for this modes then for the v\ and 1/3 stretch modes, so the primary 
excitation is seen in the umbrella mode. 

For the orientation with three hydrogens pointing towards the surface the repulsion on the 1/4 umbrella mode is 
the highest in combination with a low force constant, so this mode has a much higher excitation probability than the 
v\ and V3 stretch modes. Another interesting detail is that the 0:3 parameter is higher for the orientation with two 
than with three hydrogens pointing toward the surface, but the excitation probabilities for the ^4 umbrella in this 
orientations are equal. This is caused by a coupling between the excitation of the V2 bending and the V4 umbrella 
mode in the orientation with two hydrogens pointing towards the surface. 

We observed with Vweak (see Table [v|) that all excitation probabilities become much higher than with VMorse , except 
for the i^4 umbrella mode with two hydrogens towards the surface, which stays almost the same. It is caused by the 
fact that, although the v\ symmetrical stretch and V3 asymmetrical stretch contribute almost completely to VMorse , 
the ^4 umbrella does so just for a small part. The umbrella contributes primarily, and the V2 bending completely, 
to the harmonic terms of intramolecular PES Vi n t ra - Kveak gives only a lowering in the VM orse terms of VJntra, so 
we should expect primarily an higher excitation probability in the v\ and V3 stretch modes. This will also cause a 
higher excitation probability of the other modes, because the turn-around point will be some what later, which give 
effectively more repulsion on the other modes. 

We also observed that for both stretches the excitation probabilities shows the following trend; three hydrogens 
< two hydrogens < one hydrogen pointing towards the surface. For one hydrogen pointing towards the surface the 
excitation probability of the z/ 3 asymmetrical stretch is around twice that of the v\ symmetrical stretch, they are 
almost equal for two hydrogens, and for three hydrogens the v\ is twice the 1/3 stretch. Some of these trends can also 
be found for the PES with VMorse; but are not always that obvious. They follow the same trends as the a parameters 
of V sur{ in Table |fiL 

Kihift (see TableM) gives also higher excitation probabilities for almost all modes than VMorse, but for the v<i bending 
mode in the two hydrogens pointing towards the surface orientation it became lower. This means that the repulsion 
of the surface is here caused for a large part by the Vshift terms (see Eq. (22)), where the vi bending modes don't 
contribute. It is also in agreement with the fact that we have observed the excitation maximum earlier, so the T4 U rf 
repulsion on the V2 bending modes will be lower. Fig. ^ shows that the repulsion for the PES with Vshift is stronger 
than for the other three PESs. 

The excitations in the v\ and V3 stretch modes are extremely high with V^hift, because they contribute strongly 
to the C-H bond elongation. The trend in the excitation probabilities of the symmetrical stretch mode is caused 
by the number of bonds pointing towards the surface, because the 7^2 parameters are for all bonds and orientations 
the same, so only the strength of the coupling with the shift-factor dominates. This effect is illustrated in Fig. ||. 
We didn't observe this trend in the V3 asymmetrical stretch mode, because there is also difference between the 7^4 



parameters as can be seen in Table III and IV 



B. Structure Deformation 



If we put the methane molecule far from the surface, then all PESs are identical. For this situation we calculated a 
bond distance of 1.165 A and bond angle of 109.5 degrees in the ground state. The results of the maximum structure 



deformations are shown in the Table VI. Fig. |3| shows the names of the bonds and angles for the three orientations. 

We observed that the PESs with VMorse an d Vwcak give again the same trends, but that for the PES with Vweak 
these trends are much stronger. This is in agreement with the observations discussed for the excitation probabilities 
above. The deformations for the PES with Khift are dominated by the change of the bond distances for the bonds 
which are pointing towards the surface. These bonds become longer for all orientations in the same order. So there 



is no orientational effect and this will be probably caused completely by the Vshift terms (see Eq. (22)). 

For the orientation with one hydrogen pointing towards the surface (see Table VI), we observed that the bond 
pointing towards the surface becomes shorter, as expected, for the PES with VMorse and even more with Vwoak- This 
is caused by the repulsion of the V SUI f terms, which works in the direction of the bond axes. We observed also some 
small bond angle deformation. This is probably a secondary effect of the change in the bond distance and correlates to 
excitation of the 1/4 umbrella mode. Remarkably, for Vshift the bond length increases. Clearly the effect of the change 
in equilibrium bond length is more effective than the repulsion with the surface. The shift effect may be somewhat 
too large, but the observed change (0.232 A) is much lower than the shift of the equilibrium (0.54 A) on which we 
fitted V^hift- 

We also observed for all orientations with Vshift an earlier turn-around point, which is in agreement with the earlier 
observed excitation maximum. This is caused by the extra repulsion contribution in the V^htft terms and by the longer 
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bond length; a longer bond gives a higher repulsion for the same position of the center of mass. Fig. [j] shows this 
effect in the v% asymmetrical stretch mode. 

The orientation with two hydrogens pointing towards the surface (see Table VI) gives also shorter bond distances 
for the bonds pointing towards the surface for the PESs with VMorsc and Vweak, although very little in comparison 
to the one hydrogen orientation. We expected this, because just half of the repulsion is now in the direction of the 
bonds. The other half is perpendicular on the C-H bonds and this makes the bond angle larger. The bond angles 
between the bonds pointing towards the surface show a quite large deformation. This was already expected from the 
excitation of the vi bending and v± umbrella modes. 

We also observed that a smaller bond length of the bonds pointing towards the surface correlates with a larger 
bond angle of the angles between the bonds pointing towards the surface and a higher excitation probability of the 
v-2, bending mode. This can be explained as follows. When the bonds become shorter the center of mass can come 
closer to the surface. Because of this there will be more repulsion on the V2 bending mode, which will cause a larger 
bond angle deformation. 

Most of the deformations of the other bond angles can be explained as an indirect effect of the deformation of this 
bond angle, but the deformation of the angle between the bonds pointing away from the surface which was seen at 
the PES with Kshift must be caused by the excitation of the umbrella mode. 

We observed almost no change in bond length for the PESs with VMorsc and Vweak at the three hydrogens towards 
the surface orientation (see Table VI). Almost all energy is absorbed by the v\ umbrella mode, so this gives quite 
large bond angle deformations. The bond angle deformation is larger for the PES with Vweak than with VMorsc PES; 
in agreement with higher excitation probability for the PES with Vwcak- The PES with V^hift has around the same 
excitation probabilities for the 1^4 umbrella mode as the PES with Vweak, but the bond angle deformation is the same 
as for the PES with VMorse- This may be caused by the higher excitations of the v-i and v^, stretch modes and the 
longer bonds pointing towards the surface at the PES with V^hift, which can make it harder to deform the bond angles. 



C. Dissociation models 



Several dissociation mechanisms for direct methane dissociation on transition metals have been suggested, 
proposals can be related to two main ideas. One of them is the breaking of a single C-H bond in the initial collision! 
This noadelris most suggested in literature and also all wavepacket simulations have focussed on the effects of this 
model.EMJa The other mechanism is often called "splats" and suggests that the -critical requirement for methane 
dissociation is angular deformation of methane which allows a Ni-C bond to form.d Even though we did not try to 
describe the dissociation itself, we would like to discuss the implication of our simulation for the dissociation. 

The "splats" model seems easiest to discuss. Angular deformation is related to the V2 bending and the va umbrella 
mode. The excitation probabilities in Table [v| seem to indicate that these modes are easy to excite, and that angular 
deformation should be large. These excitation probabilities are misleading, however. Table VI shows that, although 
the excitation probabilities depend on the PES, the changes in the bond angles do not. Moreover, these changes are 
quite small. They are largest for the orientation with two hydrogen atoms pointing towards to surface, for which they 
are about 8 degrees at most. This seems much too small to enable the formation of a Ni-C bond. A previous estimate 
of the bond angle deformation gave an energy jpf 68.1kJ/mol to get three hydrogens and the carbon in one plane. 
This correspond to a bond angle of 120degrees.B We find with three hydrogens pointing towards the surface that the 
bond angle only changes to about 112 degrees at the higher energy of 96kJ/mol. The reason for this difference is that 
for the older estimate it was assumed that all translational energy can be used to deform the bond angle, whereas 
Tables [v| and VI clearly show that this is not correct. 

The excitation probabilities for the stretch modes depend strongly on the orientation of the molecule and on the 
PES. For the PESs without a change of the equilibrium bond length the excitation probabilities are only appreciably 
with one hydrogen pointing towards the surface. This is, however, not a favourable orientation for the dissociation. 
Moreover, Table [vi| shows that in those cases the repulsion with the surface shortens the C-H bond. The same 
holds for the orientation with two hydrogens pointing towards the surface, although to a lesser extent, whereas 
the bond length changes hardly at all with three hydrogens pointing towards the surface. This indicates that the 
intramolecular PES needs the bond elongation to overcome the repulsion of the surface that shortens the C-H bond, 
and to get dissociation. This agrees complet ely w ith electronic structure calculations that yield a late barrier for 
dissociation with a very elongated C-H bondr3'EJ The large excitation probabilities for the PES with the elongated 
equilibrium bond length should not be over-interpreted, however. They are to be expected even if the molecule stays 
in its (shifted) vibrational ground state. More telling is that for this PES we do find at least some inelastic scattering. 
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IV. CONCLUSIONS 



We have done wavepacket simulations on the scattering of methane from a flat Ni(lll) model surface with a fixed 
orientation with one, two, or three hydrogens pointing towards the surface. We used the MCTDH method and four 
different model PESs for each orientation. We used a translational energy of up to 96 k J/ mol and all internal vibrations 
in the ground state. The scattering was in all cases predominantly elastic. 

When the molecule hits the surface, we always observe vibrational excitations of the 1^4 umbrella and bending 
modes, especially in the orientations with two or three hydrogens pointing towards the surface. This is due to a 
favorable coupling that originates from the repulsive interaction with the surface, and the low excitation energies. 
Deformations of the molecule are predominantly in the bond angles. The changes in the bond angles are, however, 
too small to allow for the formation of a Ni-C bond, as suggested in the "splats" model of methane dissociation. 

Appreciable excitations of the v\ and ^3 stretch modes when methane hits the surface are only observed when one 
hydrogen atom points towards the surface, or when the intramolecular PES has an elongated equilibrium C-H bond 
length close to the surface. The repulsion of the surface shortens the C-H bond. This can only be overcome when the 
intramolecular PES incorporates the effect of a longer equilibrium C-H bond length caused by overlap of occupied 
surface orbitals with the antibonding orbitals of methane. This agrees with quantum chemical calculations, which 
show a late barrier for dissociation. 

The simulations with these model PESs show that the internal vibrations play an important role in the dissociation 
mechanism. Excitation probabilities when the molecule hits the surface show how the translational energy is converted 
into vibrational energy and it is distributed over the internal modes. These probabilities vary strongly with the PES. 
As only few internal vibrations contribute to the dissociation, it is important to obtain more information on the real 
PES for this system. 
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690 890 , 10p0 12p0 14p0 16p0 18p0 20pp SpO 8p0 10pp 12pp 14p0 16p0 18p0 20p0 




6p0 8p0 10p0 12p0 14p0 16p0 18p0 20p0 6p0 8p0 10pO 12p0 14p0 16p0 18p0 20p0 




FIG. 1. Contour plots in the energy range from to 0.12 Hartree of a cross-section in the orientation with one hydrogen 
atom pointing towards the surface for the translational mode X\ and the V3 asymmetrical mode X4 (an 01 mode in C3,,), 
both in atomic displacement, which show well the different behaviours of the four PESs; (a) harmonic PES, (b) an anharmonic 
intramolecular PES, (c) intramolecular PES with weakening C-H bond, and (d) intramolecular PES with elongation of the 
C-H bond. All other coordinates set to zero. 
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6qo sqo iopo 12p0 14p0 16p0 i8po 2opo 



(a) 




690 8p0 , lOpp , 12pp 14p0 16p0 18p0 20pO 



(b) 




6p0 8p0 1090 12p0 14p0 16p0 18p0 20p0 



(c) 




FIG. 2. Contour plots in the energy range from to 0.12 Hartree of a cross-section of the PES with elongation of the C-H 
bond in the translational Xi and the v\ symmetrical stretch mode X2, both in atomic displacement; (a) one, (b) two, and (c) 
three hydrogens pointing towards the surface. All other coordinates set to zero. 
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FIG. 3. Schematic representation of the i?down and R up bonds, and the adown, aside and a up angles for the three orientations: 
(a) one, (b) two, and (c) three hydrogens pointing towards the surface. 
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TABLE I. Overview of the relations between the mass-weighted coordinates Xi \ the force constants fc; (in atomic units), the 
designation, and the symmetry in Td, Cs v and Gat,. 



i 


hi 


designation 


T d 


C3d 


Ca„ 


1 




translation 


ti 


ai 


01 


2 


1.780 • 10 -4 


V\ 


symmetrical stretch 


('i 


ai 


a\ 


3 


3.599 • 10~ 5 


V4 


umbrella 


ti 


ai 


0,1 


4 


1.892 • 10 -4 


V3 


asymmetrical stretch 


ti 


ai 


ai 


5 


4.894 ■ 10 -5 


v> 


bending 


e 


e 


01 


6 


4.894 • 10 -5 


y> 


bending 


e 


e 


oa 


7 


3.599 • 10~ 5 


V4 


umbrella 


ti 


e 


bi 


8 


3.599 ■ 10~ 5 


V4 


umbrella 


ti 


e 




9 


1.892 • 10 -4 


vz 


asymmetrical stretch 


ti 


e 


bi 


10 


1.892 • 10~ 4 


vz 


asymmetrical stretch 


ti 


e 


h 



TABLE II. a and f3 values (in atomic units) of V SUI { with one, two or three hydrogens pointing towards the surface (see Eqs. 
&, & ^d (Eg)). 









one 




two 




three 


ai 


6.281 


10 _J 


6.281 


10 _J 


6.281 


10 -3 


a> 


1.256 


10~ 2 


7.252 


10 -3 


4.187 


10 -3 


«3 


4.226 


10" 3 


-7.931 


10 -3 


-1.198 


10 -2 


a 4 


-2.040 


10 -2 


-7.445 


10 -3 


-3.128 


10~ 3 


«5 






-1.026 


10~ 2 






/3i 










5.921 


10 -3 


ft 










1.026 


10~ 2 


ft 






6.079 


10 -3 


2.866 


10 -3 


& 










4.963 


10 -3 


ft 






6.476 


10" 3 


3.053 


10 -3 


ft 










5.288 


10~ 3 



TABLE III. 7 values (in atomic units) of VLiorse with one and three hydrogens pointing towards the surface (see Eq. (fig)). 





one 




three 




value 


712, 722, 732, 742 




712, 722, 732, 742 


1.079 


10 -2 


713, —3723, —3733 


-3743 


—713 , 3723 , 3733 , 3743 


1.359 


10 -3 


714, -3724, —3734 


-3744 


—714, 3724, 3734, 3744 


-1.966 


10~ 2 


717,718,719, 71,10, 


728, 72,10 


717, 718, 719, 71.10, 728, 72,10 




0.0 


727, —2737, —2747 




—727, 2737, 2747 


1.282 


10 -3 


738, —748 




738, —748 


-1.110 


10 -3 


729, —2739, —2749 




—729, 2739, 2749 


-1.853 


10 -2 


73,10, —74,10 




—73,10, 74,10 


1.605 


10 -2 



TABLE IV. 7 values (in atomic units) of Viyiorsc with two hydrogens pointing towards the surface (see Eq. (U8I)). 





two 










value 


712, 722, 732, 742 

713 , 723, —733, —743, 717, 
714, 724, —734, —744, 719, 


—727, 737, 
—729,739, 


-747, 718, - 

—749, 71,10, 


'728, —738, 748 
—72,10, —73,10, 74,10 


1.079 
-7.849 
1.135 


10-'^ 

10" 4 
10 -2 
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TABLE V. Excitation probabilities, at an initial translational energy of 96 kj/mol and all initial vibrational states in the 
ground state, for the three different PESs in the a\ modes of the C^v and C2« symmetry, with one, two or three hydrogens 
pointing towards the surface. These modes are a vi(ai) symmetrical stretch, a v%(e) bending , a ^3(^2) asymmetrical stretch, 
and a 1/4 (£2) umbrella in the Td symmetry. The PESs are: An anharmonic intramolecular PES (Morse, see Eq. (p"g|)), an 
intramolecular PES with weakening C-H bonds (weak, see Eq. (^)), and an intramolecular PES with elongation of the C-H 
bonds (shift, see Eq. (||)). 



orientation 


PES 


vi(ai) stretch 


V2(e) bending ^3 (£2) stretch 


v^(t-z) umbrella 


one 


Morse 

weak 

shift 


0.023 
0.135 
0.340 




0.054 
0.308 
0.727 


0.030 
0.075 
0.067 


two 


Morse 

weak 

shift 


0.010 
0.067 
0.707 


0.108 
0.135 
0.065 


0.009 
0.073 
0.768 


0.104 
0.107 
0.331 


three 


Morse 

weak 

shift 


0.003 
0.038 
0.819 




0.006 
0.019 
0.674 


0.102 
0.208 
0.214 


TABLE VI. Structure deformation, at an initial translational energy of 96 kj/mol and all initial vibrational states in 
the ground state, with one, two, and three hydrogen pointing towards the surface for three different PESs: An anharmonic 
intramolecular PES (Morse, see Eq. (|l8|)), an intramolecular PES with weakening C-H bonds (weak, see Eq. (pp|)), and an 
intramolecular PES with elongation of the C-H bonds (shift, see Eq. (§2])). See Fig.[] for the meaning i?down, Ru P , Qdown, a up , 
and a s idc- All bond distances are given in A and all bond angles are given in degrees. The gas phase values correspond to the 
bond distance and angles far from the surface, and are the same for all three PES types. 


Orientation 


PES 


^down 


-Rup 




gas phase 


all 


1.165 




109.5 




one 


Morse 

weak 

shift 


1.123 
1.051 
1.397 


1.166 
1.167 
1.158 


108.5 110.4 
108.2 110.7 
109.9 109.0 


two 


Morse 

weak 

shift 


1.157 
1.137 
1.386 


1.165 
1.163 
1.159 


116.7 108.9 107.8 

117.6 109.0 107.5 

115.7 106.9 108.5 


three 


Morse 

weak 

shift 


1.164 
1.167 
1.389 


1.160 
1.155 
1.144 


111.2 107.7 
112.4 106.3 
111.2 107.7 
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